Edge states and determination of pairing symmetry in superconducting Sr2Ru04 
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We calculate the energy dispersion of the surface Andreev states and their contribution to tunnel- 
ing conductance for the order parameters with horizontal and vertical lines of nodes proposed for 
superconducting Sr2Ru04. For vertical lines, we find double peaks in tunneling spectra reflecting 
the van Hove singularities in the density of surface states originating from the turning points in their 
energy dispersion. For horizontal lines, we find a single cusp-hke peak at zero bias, which agrees 
very well with the experimental data on tunneling in Sr2Ru04. 
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I. INTRODUCTION 

Numerous experiments suggest that the superconduct- 
ing state of Sr2Ru04 0| is unconventional (see review 
[ 2 1 ) . Strong suppression of Tc by non-magnetic impurities 
1 3| and absence of the Hebel-Slichter peak in NMR exper- 
iments [D indicate that it is not s-wave. The absence of 
spin susceptibility reduction below Tc in the Knight shift 
measurements indicates the spin-triplet pairing The 
/iSR experiment suggests a superconducting state with 
broken time-reversal symmetry Small-angle neutron 
scattering reveals a square vortex lattice, which is in- 
terpreted as the indication for a two-component order 
parameter These experiments initially led to sug- 
gestion of the 2D chiral (time- reversal non- invariant), 
isotropic, nodeless p-wave pairing potential How- 
ever, the power-law temperature dependences found in 
specific heat Q , nuclear relaxation rate jl^] , penetration 
depth thermal conductivity and ultrasonic at- 
tenuation (l^Jl^ indicate nodes in the energy gap. In 
response to these experiments, the following alternative 
order parameters were proposed p5[ | : anisotropic p-wave 
[ p6[ , p-wave with horizontal lines of nodes |l7j, and /- 
waves with vertical 18 1 or horizontal lines of nodes 
p9[ . It was also proposed that the a and f3 bands of 
Sr2Ru04 are either not superconducting |20| or have hor- 
izontal lines of nodes pit]. Big in- plane anisotropy of ul- 
trasound attenuation |ll4| may support the vertical lines 
of nodes. However, thermal conductivity depends very 
little on the orientation of an in-plane magnetic field , 
which is against the vertical lines |p3|| . 

Electron tunneling between a normal metal and a su- 
perconductor proved to be an important tool in determin- 
ing superconducting symmetry. Observation of the zero- 
bias conductance peak (ZBCP) |2j] due to the formation 
of the midgap Andreev bound states |2^ confirmed the 
d-wave symmetry of the high-Tc cuprates. Electron tun- 
neling into a three-dimensional p-wave superconductor 
with the pairing potential corresponding to the B-phase 
of ^He was considered in Ref. [|2^. Electron tunneling in 
Sr2Ru04 was studied theoretically in Refs. and |28] 



for the 2D isotropic chiral unitary and non-unitary p- 
waves. In this paper, we calculate the tunneling conduc- 
tance curves for the alternative order parameters listed 
above and compare them with experiments |23,M. 



II. GENERAL FORMALISM 

We model the tunneling contact by two semi-infinite 
regions, normal (N) and superconducting (S), with a flat 
interface (I) perpendicular to the a axis of Sr2Ru04 (see 
Fig. |l|). The x and y axes are selected along the crys- 
tal axes a and 6, respectively. The tunneling conduc- 
tance is calculated at zero temperature by solving the 
Bogoliubov-de Gennes equations in the ballistic regime 
following Refs. [ ^Ip^ . The tunneling barrier is mod- 
eled by a delta- function potential of strength H, so the 
boundary conditions at the interface are ipn\i = tps\i and 
Vn4'n\i = Vgipsli — '^i'H'ipnli- Here -ipn and ips are the elec- 
tron wave functions in the normal and superconducting 
regions, and Vn,s are the velocity components perpendic- 
ular to the interface. The pairing potential of a triplet su- 
perconductor can be expressed as A = ioyicr ■ d)A(k, x), 
where <t are the Pauli matrices operating on the electron 
spin indices, the vector d indicates the direction of spin 
polarization, and k is the relative momentum of elec- 
trons in a Cooper pair. In this paper, we consider the 
cases where d does not depend on k and is pinned to the 
c axis of Sr2Ru04. The spatial dependence of the pairing 
potential is taken to be step-like: A(k, a;) — A(k) 8(a;). 
Assuming that the electron momentum ky parallel to the 
interface is conserved, the tunneling conductance per one 
layer, G(y), can be written as an integral over ky |^ , ^ : 



Giy) 

GiV,ky) = D 



i + |r|2D-(i-i:))|r|4 



1 - (1 - Z?) r2 exp(i$) 



(1) 

(2) 



In Eq. ([^), e is the electron charge, h is the Planck 
constant, I is the length of the interface, and V is the 
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FIG. 1. Fermi surfaces of the normal (N) and supercon- 
ducting (S) materials forming the interface I. The points A 
and B on the Fermi surface S are connected by specular re- 
flection from the interface. 



bias voltage. We consider the case where the Fermi 
surfaces of both the normal metal and superconductor 



are circular, and their radii are fc^^ > kp' , as shown 
in Fig. |l|. The limits of integration are set by the 



(s) 

smaller Fermi momentum: k^^^ = kp . In Eq. ^ 
D{ky) = 4,Vn{ky)vs{ky)/{[vn{ky) + v,(ky)f + is the 
normal-state transmission coefficient, where Vn.s{ky) are 
the Fermi velocities components perpendicular to the in- 
terface, ^{ky) ~ (j>A{ky) — 4>B{ky) is the phase difference 
of the superconducting pairing potentials at the points on 
the Fermi surface connected by electron reflection from 
the interface (points A and B in Fig. |l|) [Q. T{V, ky) is 



eV - sgn(el/)^(el/)2- |A(fcj,)| 



|A(fc, 



r = 



zV - i^\^{ky)\^ - {evy 
|A(fc,)| 



\eV\ > |A|, 
eV\ < |A|. (3) 



Denoting G{V,ky) as Q{V,ky) for the subgap voltage 
\eV\ < |A|, where |r(y)| = 1, we rewrite Eq. (||) as 



GiV,ky) 



DHky)/2[l-Diky)] 
D^{ky)/4[l-D{ky)]+F^{V,ky)' 



(4) 



F = y/l - |ey/A|2cos($/2) - (ey/|A|) sin($/2). (5) 

For a high barrier, D sa VnVs/H'^ <C 1. Then Eq. (^) 
becomes g{V, ky) w TTD{ky) S[F{V, ky)], and Eq. (0) gives 



giv)^c 



/ '-^S[e-Eiky)]^j: 



)i^(#)r 



(8) 



where fc^"*-* is the j-th root of the equation E{ky) = e. 

The delta-function in Eq. (^) indicates that the subgap 
tunneling takes place when the bias voltage matches the 
energy of a surface state: eV — E{ky). Denoting the 



j-th root of this equation as ky 
function in Eq. (^, we find 



0) 



and resolving the delta- 



Q{v)^cY^ 



)V-{evy^ 
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(9) 



Because the denominators in Eqs. (||) and (||) are the 
same, both DOS and tunneling conductance exhibit 
peaks (the inverse-square-root van Hove singularities) 
at the turning points of the energy dispersion, where 
dkyE{ky) = 0. The peak positions depend solely on 
the pairing potential A(fcy) through Eq. (|^) and not on 
the band structure details. The munbcr of such turning 
points does not change upon small continuous deforma- 
tion of A (fcy), so it is a topological feature. However, 
because of the additional factors in Eq. (^, tunneling 
conductance is not simply proportional to DOS, as of- 
ten assumed. Eq. (^) vanishes at ky kp\ where the 
normal component Vg of the Fermi velocity goes to zero. 
It also vanishes when A{ky) ^ or E(ky) ±A{ky) 
(i.e. ^{ky) 0, 2tt). In these cases, the quasipar- 
ticle localization length perpendicular to the interface, 
A = hvs/ ^|Ap — E"^, diverges |Q; thus, the probability 
to find a quasiparticle on the surface vanishes. 



III. APPLICATION TO SraRuO* 

In Figs. H ^, ^, and ||, we show how this general for- 
malism applies to the alternative superconducting order 
parameters proposed for Sr2Ru04 in the literature. In 
order to focus on the features of superconductor, rather 



y 



than normal metal, we took k 



(n) 
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5k'p'. Then, the per- 
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dky Vn{ky) Vs{ky) 5[F [V , ky)], (6) pendicular velocity of the normal metal is approximately 



constant for all momenta |fc„ 



< 4^^ 



,(") 



while 



where C = 2TTe l/h 7i . The delta-function contributes 
to the integral d) when F{V, fc^) = 0. Using Eq. (|), this 
condition can be written as eV = E{ky), where 



£;(fc,) = |A(/c,)|cos[<i>(fc,)/2] 



(7) 



for sin($/2) > 0, i.e. < $ < 2tt, which is the appropri- 
ate interval for <I>. E{ky) (|^) is nothing but the energy of 
a surface Andreev state with the momentum ky obtained 
for the impenetrable barrier (Ti. — > oo) p5[ |. The energy 
density of these surface states (DOS) is 



in the superconductor it is Vs — Vp'^ 1 — [ky/k^p^)'^, 

where and v^p' are the corresponding Fermi veloc- 
ities. The numerical calculations are done for the high 

barrier Z = Ti./ \l v'^p\'^p' = 3. As the figures illustrate, 
in this limit, the approximate curves given by Eq. (|^), 
which includes only the contribution of the localized sur- 
face states, agree well with the exact curves obtained 
from Eqs. (|l]) and (||), which also include the contribu- 
tion of the extended bulk states. 
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FIG. 2. Left panel: Energy dispersion E\(ky) mm for the 
2D isotropic p-wave. Right paneh The corresponding exact 
[Eqs. (0) and (^), dash-dotted hne] and approximate [Eq. (p^, 
dotted hne] tunnehng conductance curves; DOS (sohd hne) 
and tunnehng conductance (long-dashed line) for the p-wave 
with horizontal lines of nodes [Eq. (p^j. 



A. Nodeless pairing potential 

For the 2D isotropic chiral p-wave, the pairing poten- 
tial and the energy dispersion of the edge states are [p|j28| 



Al(k) = A0(fc. + iky)lkf, El(ky) = A0ky/k^;\ (lO) 

Substituting Eq. (|l^) into Eqs. (|) and (|), we find 



k^F^/hAo 



2Go[l 



= Po, 
[eV/Aor] 



(11) 
(12) 



where Go = Tre^lkp^v^p\p^ /h'^Tl'^ is the normal-state 
conductance. The energy dispersion Ei{ky) (nffl) has no 
turning points, and the corresponding DOS ( |ll[ ) is flat: 
p(e) = const. The tunneling conductance curve G{V) 
( |l^ ) is parabolic for \eV\ < Aq, attaining the maximal 
value 2Go at V" = and vanishing at \eV\ — Aq, as 
shown in Fig. ^. For \eV\ > Aq, G{V) increases with in- 
creasing \V\ and saturates at the normal-state value Gq 
at \eV\ ^ Aq. This behavior is in agreement with pre- 
vious numerical calculations ||2^ , |2^ . One can recognize 
that the curves shown in Fig. 4 of Ref. and the cor- 
responding figures of Ref. |2^ for \eV\ < Aq are actually 
distorted parabolas. The deviations are due to the finite 
values of Z and to the finite size of the acceptance cone 
utilized in the numerical calculations |2^j2^. Our ana- 
lytical expression (|lj) applies in the case Z 3> 1 and for 
the full acceptance cone. 



B. Pairing potentials with vertical lines of nodes 

In this subsection, we consider pairing potentials with 
vertical lines of nodes, which originate from modulation 
of the gap function A(k) in the {k^, ky) plane. 

For the anisotropic chiral p-wave proposed in Ref. [p^ , 
the pairing potential and the energy dispersion are 



A2(k) = Ao[sin(fc^a/?i) 
E2{ky) = Aq sin(fcya/?i), 



i sin{kya / h)], 



(13) 
(14) 
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FIG. 3. Left panel: Energy dispersion E2{ky) [Eq. (Q, 
solid line] and gap \A2{ky)\ [Eq. (P), dashed hne] for the 
anisotropic p-wave. Right panel: the corresponding DOS 
[dashed-dotted line] and the exact (solid line) and approxi- 
mate [dashed line] tunneling conductance curves. 



where a is the lattice constant, and k 



\ ' ,-7 



0.97rn. Sub- 



stituting Eqs. (|T|) and (|T|) into Eqs. (§), (|), (|), (|) 
and (|8|), we calculate the tunneling conductance and 
DOS curves shown in Fig. ^. They exhibit peaks at 
eV ~ ±Ao originating from the turning points in the 
edge states dispersion. Equation eV = E{ky) has one 
solution ky"* for \eV\ < 0.31 Aq and two solutions for 
0.31 Aq < |eT^| < Aq. The switch causes the discontinu- 
ity of DOS and the slope change of tunneling conductance 
at eV = ±0.31 Aq. 

For the /-wave proposed in Ref. pTt], we have 



A3(k) = Ao {kl 
Es{ky) ^Ao \kl 



^1) 



[k 

k^, I ky I {k^p ^ ) . 



(15) 
(16) 



As shown in Fig. ^, the edge states dispersion has cusps at 
ky — ±fc^^/-\/2, where A3(k) vanishes, and two turning 
points in between. Consequently, the tunneling conduc- 
tance and DOS curves have peaks at eV = ±0.26 Aq. 
For another /-wave proposed in Refs. |p^,na], we have 



A4(k) 

Eiiky) 



2Ao hj,ky {k. 



iky) I (fc^'*)'^, 



2Ao ki ky/{k^;^ f. 



(17) 
(18) 



Eqs. (|T^ and ([T|) transform into Eqs. (|5|) and (|l|) upon 
a 7r/4 rotation in the x-y plane. As shown in Fig. the 
turning points in the energy dispersion and the peaks 
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FIG. 4. The same curves as in Fig. y| for the /-wave pairing 
potential (^5|) with vertical lines of nodes. 
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FIG. 5. The same curves as in Fig. |3| for the /-wave pairing 
potential (|l7[) with vertical lines of nodes. 



in the tunneling conductance and DOS curves occur at 
eV — ±0.76 Aq. Their position is approximately the 
same as in Fig. ^ for the anisotropic p-wave. However, in 
Fig. 1^, unlike in all other figures, tunneling conductance 
vanishes at zero bias, even though DOS remains finite. 
This characteristic feature of the /-wave potential ( p^ ) 
is a consequence of the square- root factor in Eq. (^) . 

Thus, tunneling conductance G(y) exhibits two peaks 
located symmetrically around y = for all pairing po- 
tentials with vertical lines of nodes described in this sub- 
section. This conclusion is in agreement with recent nu- 
merical calculations IpSl. 



C. Pairing potentials with horizontal lines of nodes 

Thus far we considered pairing potentials that do not 
depend on fcz, the momentum component perpendicular 
to the Sr2Ru04 planes. Ref. also proposed the p- 
wave with horizontal lines of nodes at certain fc^. Its 
pairing potential can be obtained from Eq. ( p^ by the 
following substitution 



Ao Aq{p^) = Aocos(p2), 



(19) 



where we switched to the dimensionless variable p 
kzc/h (c is the interplane distance). 

Now DOS is obtained by averaging p(e) from Eq. (|l 
over pz using Eq. (|l^) : 



dpz 



27r /i|Ao(p^ 



|Ao(p.)|>|e| 



(20) 



(21) 



In Eq. (|2C|), IS.q{pz) vanishes linearly withp^ at the nodes, 
which results in a logarithmic divergence of p{e) at e — > 0: 



/9(e) cx ln|Ao/e|, for |e| < Ao 



(22) 



The logarithmic divergence is cut off at |Ao(p2)| = |e|, 
because Eq. (^l]) applies only to the subgap states with 
|e| < \Kq{pz)\. The plot of p{t) is shown in Fig. H 



Tunneling conductance Q (V) is obtained by substitut- 
ing Eq. (noh into Eq. dl2) and averaging over pz : 



g{V) = 2Go 



|Ao(p.)|>|eV| 
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(23) 



(24) 



In Eq. (|23|), the divergence of J_/|Ao(pz)| at the node is 
cut off at \eV\, because Eq. ( |l2|) apphes only for \eV\ < 
|Ao(pz)|. To treat the divergence, it is convenient to 
change the variable of integration from pz to A: 



l^^l dA 
A' 



4|ey|^ 
2n 



Ao 



\eV\ 



dA 
A^ 



(25) 



Here the factors of 4 come from the two nodal lines and 
the integration over the two sides of each line. The slope 
A' = \dAo{pz) / dpz\ is approximately constant at the 
nodes: A' « Aq. Then the integrals in Eq. ( p5|) can 
be taken, and we find 



g{V) « 2Go 1 - - 



eV 
A'o 



for \eV\ < Ao. (26) 



Thus, tunneling conductance has a cusp at zero bias. The 
curve for G{V) calculated exactly starting from Eqs. (||) 
and (H) is shown in Fig. || and indeed demonstrates the 
cusp in agreement with Eq. (|26|). 

Notice the difference in Figrgbetween the singular log- 
arithmical divergence of DOS and the triangular-shaped 
cusp in tunneling conductance at zero bias. Both features 
result from the pile-up of surface Andreev states at zero 
energy caused by vanishing gap at the nodes. These are 
robust features of the chiral pairing potentials with hor- 
izontal lines of nodes, independent of the band structure 
details. While Eqs. ( ^l] ) and ( p^ ) are derived specifically 
for the gap ( p^ b eing a cosine function, the asymptotic 
Eqs. (22) and(M) are valid for any generic chiral pairing 
potential with horizontal lines of nodes. For example, the 
/-wave with horizontal lines of nodes p9| , 

A6(k) = lAo sinipz) (fc. + ikyf I (kff, (27) 
E^(ky) = sgn(fcj,) Ao \2(kylkf f - 1] I sin(p,)|, (28) 

produces similar curves for the in-plane tunneling. How- 
ever, being an odd function oipz^ it also exhibits a ZBCP 
in the c-axis tunneling. 



IV. COMPARISON WITH EXPERIMENT 

Point-contact tunneling spectroscopy |^9| found a 
ZBCP in Sr2Ru04. The data were fitted using the 2D 
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FIG. 6. Long-dashed line: Theoretical tunneling conduc- 
tance curve from our Fig. ^ calculated for the p-wave with 
horizontal lines of nodes [Eq. (jl^)] using Aq = 0.8 meV. Solid 
line: Experimental curve from Fig. 2 of Ref. jsoj for 0.8 K. 

isotropic p-wave (0) and assuming a narrow acceptance 
cone, i.e. a small value of fc™'*'^. In this situation, only the 
states with ky Ri and hence E{ky) « are probed; thus 
all pairing potentials, except the /-wave (^|), would show 
a ZBCP. However, such a ZBCP, in contrast to nonchiral 
cuprates and organic superconductors jsij , does not 
identify the pairing symmetry uniquely. 

The more recent tunneling experiments |30| ] were per- 
formed on cleaved junctions of Sr2Ru04 with inclusions 
of Ru. Although nominally their setup corresponds to 
tunneling along the c axis, the experimentalists believe 
that the actual tunneling takes place in the (a, b) plane 
via the inclusions of Ru. They found a cusp-like peak at 
zero bias in the 1.4-K phase and a sharper ZBCP in the 
3-K phase. The latter was attributed to a nonchiral su- 
perconductivity developing at the grain boundaries [ p5[ . 
An experimental curve from Ref. | |30[ | representing the 
1.4-K phase is shown in Fig. |^. It has the characteristic 
triangular shape with a cusp at zero bias and agrees very 
well with our theoretical curve calculated for the pairing 
potential ([l9|). To making the comparison, we treated 
the horizontal and vertical scales and the vertical offset 
(the background) of the theoretical curve as fitting pa- 
rameters. The best fit corresponds to the gap Aq ~ 0.8 
meV. 

The good agreement between the theoretical and ex- 
perimental curves is a strong indication in favor of the p- 
or /-wave pairing potential with horizontal lines of nodes 
in the 1.4-K phase. Further measurements showing the 
absence or presence of a ZBCP in the c-axis tunneling 
could discriminate between the p- and f-wave cases. Ex- 
periment |3^] found that an additional sharp ZBCP with 
the width of the order of 0.1 meV develops on top of the 
cusp at lower temperatures 0.5 and 0.32 K. It may be 
attributed to a contribution from the c-axis tunneling, if 
the pairing potential is the /-wave (p7|). 

The described above scenario corresponds to the hori- 
zontal nodes in the main electron band 7 of Sr2Ru04. In 
the alternative scenario proposed in Ref. ||2l[], the main 
band 7 is nodeless, but the two other bands a and P have 



horizontal lines of nodes. In this case, one would expect a 
superposition of a parabolic curve for the nodeless band 
and triangular curves for the bands with nodes. Perhaps 
the sharp ZBCP developing at low temperatures 0.5 and 
0.32 K could be interpreted as the triangular curve cor- 
responding to the gap of the order of 0.1 meV in the a 
and f3 bands. 

V. CONCLUSIONS 

We derived the analytical formula (||) for subgap tun- 
neling at low transparency of the barrier, which takes 
into account only the contribution of the surface An- 
dreev states. This formula produces tunneling curves in 
a simple and physically transparent way. For the chiral 
pairing potentials with vertical lines of nodes, tunnel- 
ing curves show double peaks, which originate from the 
turning points in the energy dispersion of the surface An- 
dreev states. On the other hand, for the chiral pairing 
potentials with horizontal lines of nodes, we find a single 
triangular-shaped peak with a cusp at zero bias, which 
results from the pile-up of the surface Andreev state at 
zero energy caused by vanishing gap at the nodes. 

Double peaks in tunneling conductance were observed 
experimentally in some point contacts with Sr2Ru04 in 
Ref. However, such double-peaked spectra were as- 
sociated with contacts with high transparency of the bar- 
rier ]29| ], a regime that is not addressed in our paper. 
On the other hand, no double peaks were observed in 
experiment pO| , which only found a single triangular- 
shaped peak at zero bias. As shown in Fig. ^, this peak 
is very well fitted by our calculations for horizontal lines 
of nodes. It is also visually similar to the single peak 
found for contacts with low transparency in Fig. 3 of ex- 
periment (which shows resistance, rather than con- 
ductance). Thus, we conclude that the superconducting 
pairing potential in Sr2Ru04 most likely has horizontal 
lines of nodes. 

Note added in proof. Quantitative comparison shows 
that our calculation for horizontal lines of nodes also fits 
the experimental data of Ref. 29 very well. The curves 
are shown at 

http://www2.physics.umd.edu/~yakovenk/talks/Sr2Ru04/ 

We are grateful to Ying Liu and Gemot GoU for send- 
ing us the experimental data of Refs. 30 and 29, and 
to Yoshiteru Maeno, Ying Liu, and Igor Zutic for useful 
discussions. This work was supported by the Packard 
Foundation and by the NSF Grant DMR-98 15094. 



5 



[1] Y. Maeno, H. Hashimoto, K. Yoshida, S. Nishizaki, T. 
Fujita, J. G. Bednorz, and F. Lichtenberg, Nature (Lon- 
don) 372, 532 (1994). 

[2] Y. Maeno, T. M. Rice, and M. Sigrist, Physics Today 54, 
#1, 42 (January 2001); 54, #3, 104 (March 2001). 

[3] A. P. Mackenzie, R. K. W. Haselwimmer, A. W. Tyler, 

G. G. Lonzarich, Y. Mori, S. Nishizaki, and Y. Maeno, 
Phys. Rev. Lett. 80, 161 (1998); Z. Q. Mao, Y. Mori, and 
Y. Maeno, Phys. Rev. B 60, 610 (1999). 

[4] K. Ishida, Y. Kitaoka, K. Asayama, S. Ikeda, S. 

Nishizaki, Y. Maeno, K. Yoshida, and T. Fujita, Phys. 

Rev. B 56, R505 (1997). 
[5] K. Ishida, H. Mukuda, Y. Kitaoka, K. Asayama, Z. Q. 

Mao, Y. Mori, and Y. Maeno, Nature (London) 396, 658 

(1998) ; K. Ishida, H. Mukuda, Y. Kitaoka, Z. Q. Mao, H. 
Fukazawa, and Y. Maeno, Phys. Rev. B 63, 060507(R) 
(2001). 

[6] G. M. Luke, Y. Fudamoto, K. M. Kojima, M. I. Larkin, J. 

Merrin, B. Nachumi, Y. J. Uemura, Y. Maeno, Z. Q. Mao, 

Y. Mori, H. Nakamura, and M. Sigrist, Nature (London) 

394, 558 (1998). 
[7] T. M. Riseman, P. G. Kealey, E. M. Forgan, A. P. 

Mackenzie, L. M. Galvin, A. W. Tyler, S. L. Lee, C. 

Ager, D. McK. Paul, C. M. Aegerter, R. Cubitt, Z. Q. 

Mao, T. Akima, and Y. Maeno, Nature (London) 396, 

242 (1998); P. G. Kealey, T. M. Riseman, E. M. Forgan, 

L. M. Galvin, A. P. Mackenzie, S. L. Lee, D. McK. Paul, 

R. Cubitt, D. F. Agterberg, R. Heeb, Z. Q. Mao, and Y. 

Maeno, Phys. Rev. Lett. 84, 6094 (2000). 
[8] T. M. Rice and M. Sigrist, J. Phys. Condens. Matter 7, 

L643 (1995); M. Sigrist, D. Agterberg, A. Furusaki, C. 

Honerkamp, K. K. Ng, T. M. Rice, and M. E. Zhito- 

mirsky, Physica C 317-318, 134 (1999). 
[9] S. NishiZaki, Y. Maeno, and Z. Mao, J. Phys. Soc. Jpn. 

69, 572 (2000). 

[10] K. Ishida, H. Mukuda, Y. Kitaoka, Z. Q. Mao, Y. Mori, 

and Y. Maeno, Phys. Rev. Lett. 84, 5387 (2000). 
[11] I. Bonalde, B. D. Yanoff, M. B. Salamon, D. J. Van Har- 

lingen, E. M. E. Chia, Z. Q. Mao, and Y. Maeno, Phys. 

Rev. Lett. 85, 4775 (2000). 
[12] M. A. Tanatar, S. Nagai, Z. Q. Mao, Y. Maeno, and T. 

Ishiguro, Physica C 341-348, 1841 (2000); Phys. Rev. B 

63, 064505 (2001). 
[13] H. Matsui, Y. Yoshida, A. Mukai, R. Settai, Y. Onuki, 

H. Takei, N. Kimura, H. Aoki, and N. Toyota, Phys. Rev. 
B 63, 060505(R) (2001). 

[14] C. Lupien, W. A. MacFarlane, C. Proust, L. Taillefer, Z. 
Q. Mao, and Y. Maeno, Phys. Rev. Lett. 86, 5986 (2001). 

[15] Vertical and horizontal lines of nodes were discussed for 
cuprate superconductors by R. J. Radtke, A. I. Liechten- 
stein, V. M. Yakovenko, and S. Das Sarma, Phys. Rev. B 
53, 5137 (1996). 

[16] K. Miyake and O. Narikiyo, Phys. Rev. Lett. 83, 1423 

(1999) . 

[17] Y. Hasegawa, K. Machida, and M. Ozaki, J. Phys. Soc. 

Jpn. 69, 336 (2000). 
[18] M. J. Graf and A. V. Balatsky, Phys. Rev. B 62, 9697 

(2000) . 

[19] H. Won and K. Maki, Europhys. Lett. 52, 427 (2000). 
[20] D. F. Agterberg, T. M. Rice, and M. Sigrist, Phys. Rev. 
Lett. 78, 3374 (1997). 



[21] M. E. Zhitomirsky and T. M. Rice, Phys. Rev. Lett. 87, 
057001 (2001). 

[22] M. A. Tanatar, M. Suzuki, S. Nagai, Z. Q. Mao, Y. 
Maeno, and T. Ishiguro, Phys. Rev. Lett. 86, 2649 
(2001); K. Izawa, H. Takahashi, H. Yamaguchi, Y. Mat- 
suda, M. Suzuki, T. Sasaki, T. Fukase, Y. Yoshida, R. 
Settai, and Y. Onuki, ibid. 86, 2 653 (2001). 

[23] T. Dahm, H. Won, and K. Maki, |cond-mat/0006301 

[24] S. Kashiwaya, Y. Tanaka, M. Koyanagi, H. Takashima, 
and K. Kajimura, Phys. Rev. B 51, 1350 (1995); M. Cov- 
ington, M. Aprili, E. Paraoanu, L. H. Greene, F. Xu, J. 
Zhu, and C. A. Mirkin, Phys. Rev. Lett. 79, 277 (1997); 
J. Y. T. Wei, N.-C. Yeh, D. F. Garrigus, and M. Strasik, 
ibid. 81, 2542 (1998). 

[25] C.-R. Hu, Phys. Rev. Lett. 72, 1526 (1994); J. Yang and 
C.-R. Hu, Phys. Rev. B 50, 16766 (1994). 

[26] L. J. Buchhohz and G. Zwicknagl, Phys. Rev. B 23, 5788 
(1981). 

[27] M. Yamashiro, Y. Tanaka, and S. Kashiwaya, Phys. 
Rev. B 56, 7847 (1997); M. Yamashiro, Y. Tanaka, Y. 
Tanuma, and S. Kashiwaya, J. Phys. Soc. Jpn. 67, 3224 
(1998); N. Yoshida, Y. Tanaka, J. Inoue, and S. Kashi- 
waya, ibtd. 68, 1071 (1999); M. Yamashiro, Y. Tanaka, 
N. Yoshida, and S. Kashiwaya, ibid. 68, 2019 (1999). 

[28] C. Honerkamp and M. Sigrist, J. Low T. Phys. Ill, 895 
(1998). 

[29] F. Laube, G. GoU, H. v. Lohneysen, M. Fogelstrom, and 

F. Lichtenberg, Phys. Rev. Lett. 84, 1595 (2000). 
[30] Z. Q. Mao, K. D. Nelson, R. Jin, Y. Liu, and Y. Maeno, 

Phys. Rev. Lett. 87, 037003 (2001). 
[31] G. E. Blonder, M. Tinkham, and T. M. Klapwijk, Phys. 

Rev. B 25, 4515 (1982). 
[32] Y. Tanaka and S. Kashiwaya, Phys. Rev. Lett. 74, 3451 

(1995). 

[33] N. Stefanakis, J. Ph ys. Condens. Matter 13, 3643 (2001); 

|cond-mat/0109498 . 
[34] K. Sengupta, I. Zutic, H.-J. Kwon, V. M. Yakovenko, and 

S. Das Sarma, Phys. Rev. B 63, 144531 (2001). 
[35] M. Sigrist and H. Monien, J. Phys. Soc. Jpn. 70, 2409 

(2001). 



6 



